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An overdamped dynamical system, biased by an external constant force, does not exhibit negative 
mobility. However, when the system is coupled to its copy, negative mobility can arise. We show it 
by the example of an experimentally realizable system of two coupled resistively shunted Josephson 
junctions. The first junction is dc-biased by a constant current and ac-driven by an unbiased 
harmonic current. In the proposed setup one can control the transport properties of the second 
junction by manipulating the first junction only. We demonstrate that the second junction can 
exhibit both absolute negative resistance (ANR) near zero bias and negative resistance in the non- 
linear regimes (NNR). These anomalous transport features occur for restricted windows of the 
coupling constant. 



Systems under non-equilibrium conditions can display 
new features as well as unexpected phenomena and pro- 
cesses which are forbidden in equilibrium systems by fun- 
damental laws. One of the most prominent examples in- 
clude the phenomena of negative mobility (conductance, 
resistance) [lj: when a constant force is applied to a mo- 
bile particle, it moves in the direction opposite to that 
of the force. This phenomenon is impossible in equilib- 
rium states because it would violate the second law of 
thermodynamics. A simple model is formulated in terms 
of the Newton equation describing the one-dimensional 
motion of a classical Brownian particle. In the minimal 
model, the particle moves in a symmetric spatially pe- 
riodic potential, is biased by a static force, and driven 
by the unbiased harmonic force ^4cos(co>£) [2H4]. This 
system is out of equilibrium and displays both absolute 
negative mobility (ANM) around zero static applied force 
(the linear response regime) and negative mobility in the 
non-linear response regime (NNM) . In the linear response 
regime, the long-time averaged particle velocity tends to 
zero when the static force tends to zero. In the nonlinear 
response regime, the particle velocity can tend to zero 
even if the static force is far from zero. It is known that 
the corresponding overdamped dynamical system does 
not exhibit negative mobility: the inertial term in the 
Newton equation is absolutely necessary for the negative 
mobility to occur. However, phenomena which are ab- 
sent in a single element perhaps can occur in a system 
of coupled elements. The physics of many-body systems 
provides hundreds of examples. Therefore, we are going 
to check whether it is possible to construct a system of 
two coupled overdamped particles that exhibits negative 
mobility. We can speculate that the combined effect of 
the time-dependent driving and the interaction between 
the two particles can radically modify a one-particle over- 
damped dynamics yielding transport anomalies like ANM 
or NNM. We apply the static and time-periodic forces to 
only one (say the first) particle and observe the response 
of the second particle. In a more general context, one 
can manipulate the first particle and by this means con- 



trol the transport properties of the second particle. As 
a possible realization of such a system we exploit two 
Josephson junctions which play the role of the two par- 
ticles in the mechanical framework. 

Model -In a more general case, we can consider a sys- 
tem consisting of two coupled resistively shunted Joseph- 
son junctions with the critical currents (I c i ,I C 2), resis- 
tances (Ri,R 2 ) and phases (^1,^2), respectively. The 
system is externally shunted by the resistance R3 as 
shown in Fig. [T] Moreover, the junctions can be dc- 
biased by the currents I\ and I 2l and ac-driven by the 
currents a\ cos(u;£) and a 2 cos(ut), respectively. We con- 
sider the small junction area limit and a regime where 
photon-assisted tunneling phenomena do not contribute. 
The phase dynamics of the junctions is determined by 
two equations, the dimensionless form of which read [5 J 

0i = ii — I c i sin <pi + bi cos(ujt) 

+ ol (J 2 - I c2 sin <t> 2 ) + VD £'(*), (la) 
(j) 2 = a /3 (I 2 - I C 2 sin (j) 2 ) + b 2 cos(ut) 

+ a(h- I cl sin 0x) + Vc^D £"(t), (lb) 




FIG. 1: The system of two resistively shunted Josephson junc- 
tions coupled by an external shunt resistance R3 and driven by 
the currents h(t) — h+ai cos(u;£) and h(t) = I2+CL2 cos(ujt). 
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where the dot denotes a time derivative, and the param- 
eters (3 = I + R3/R1, bi = di +aa2 and 62 = a(ai +^2). 
The coupling parameter a = ^2/(^2 + R3) £ [0, 1] can 
be changed by adjusting the external resistance R%. We 
use the same dimensionless units as Nerenberg et al. [5]: 
the current unit is the averaged critical supercurrent 
Ic = {Id + ^c2)/2, and the time unit is fi/2eVb, where 
Vb = I c i^i (i?2 + i?3) /{Ri + R2 + R3) is the characteris- 
tic voltage. We assume that all resistors are at the same 
temperature T, and that the noise sources (the John- 
son thermal noise) are represented by zero-mean white 
noises £i{t) {i = 1,2,3) which are delta-correlated, i.e. 
(€i{i)£j{s)) = S i:j S{t - s) for ij e {1,2,3}. The noises 
£'(£) and £"(£), which appear in Eqs. ([I]), are linear com- 
binations of the noises = 1,2,3. The dimension- 
less noise strength is D = AeksT /hl c . Eqs. are the 
extended version of the system studied by Nerenberg et 
al. [SJ [6] which now include ac driving and noise terms. 

Eqs. ([I]) describe over damped dynamics of a hypothet- 
ical mechanical system of two interacting particles of 
coordinates x\ = <j)\ and X2 = 4>2, respectively. The 
phase space of the deterministic system ([I]) is three- 
dimensional, namely {x\ = 4>i,X2 = </>2,^3 = wt}. 
Three is the minimal dimension of the phase space to dis- 
play chaotic evolution which is an important feature for 
anomalous transport to occur [2H3]. At a non-zero tem- 
perature, D > 0, the Johnson thermal fluctuations acti- 
vate a diffusive dynamics where stochastic escape events 
among existing attract or s become possible. Moreover, 
the system can now visit any part of the phase space and 
evolve within some finite time interval by closely follow- 
ing any existing orbits, either stable or unstable. 

Driving applied to one junction. - The detailed anal- 
ysis of the full system ([I]) will be presented elsewhere. 
Here, we consider a simplified system of two identical 
junctions {Ri = i?2, Id = Ic2 = 1, aft = 1) and the 
case when I 2 = and 02 = - that is, when only the 
first junction is dc- biased and ac-driven. In this way, 
we want to manipulate the first junction only and ob- 
serve the response of the second junction. We ask under 
what conditions we can control transport properties of 
the second junction. To this aim, we numerically study 
the dimensionless long-time averaged voltage V2 = (^2) 
across the second junction. The corresponding asymp- 
totic averaged voltage across the first junction is denoted 
as vi = (0i ). The long-time physical voltage is then ex- 
pressed as Vi — VoVi {i = 1,2). If the problem is formu- 
lated in terms of overdamped motion of classical Brown- 
ian particles, the voltage Vi translates into the averaged 
velocity of the first or second particles, respectively, and 
the dc current I\ translates into an external constant 
force acting only on the first particle. The junction resis- 
tance (or equivalently conductance) translates then into 
the particle mobility. One can use this analogy to sim- 
plify the visualization of transport processes in junctions. 
The voltage V{ = Vi{I\) is typically a nonlinear and non- 



monotonic function of the dc current I\. In the "normal" 
transport regime, the voltage is positive Vi > for pos- 
itive bias I\ > 0, i.e. the "nonlinear resistance" or the 
static resistance R = R{Ii) — Vi{I\)jI\ at a fixed bias 
current 1\ is positive. If R < 0, i.e. when the system 
response is opposite to the external load, it refers to the 
anomalous transport regime with ANR or NNR. 

We begin the analysis of the system by some gen- 
eral remarks about its long-time behavior. As expected 
from symmetry, in the zero bias case, i.e. when I\ = 0, 
there is no net transport. If the dc bias I\ is sufficiently 
large in comparison to the amplitude a\ of the ac driving, 
the voltage across both junctions has the same sign as the 
dc current. This is rather obvious because the driving is 
not important in such a regime. More interesting effects 
can take place in the regime of small 1\. However, the pa- 
rameter space {a, Ji, ai, a;, D} is 5-dimensional and thus 
too large for an extensive numerical scan. Therefore, 
a number of low resolution scans over the parameters 
{a : ai : u} at fixed values of I\ and D = was first per- 
formed. Next, the interesting regions of the parameter 
space were analyzed in more detail and at a higher res- 
olution. These initial scans were done for initial values 
of the phases <pi and (j)2 randomly chosen from the in- 
terval [0; 2tt]. They were then supplemented by scans at 
non-zero temperatures. We have found that for Ii < 0.1 
the modulus of the average voltage across both junctions 
takes its highest values for uj G (0, 1) and is gradually 
diminished for higher frequencies. At large ac-driving 
frequencies {u > 5) there is no noticeable net transport 
regardless of the values of other parameters. 

For fixed but small values of the dc bias Ii < 0.1, strips 
of non-zero average voltage are clearly visible in the pa- 
rameter space {ai,co>, a}. For weak coupling of two junc- 
tions (small a), there is no net transport in the second 
junction. For stronger coupling, strips of non-zero aver- 
age voltage start to appear at progressively lower values 
of the amplitude a\ of the ac driving. The strips are also 
visible in the plots of the average voltage of the first junc- 
tion, which means that they represent regimes of the pa- 
rameter space where both junctions operate in synchrony. 
E.g. if h = 0.05, for a e [0.66; 0.9] and a x ~ 2.0 ± 0.5 
the structures visible in the plots become more complex 
and regions of negative conductance of the second junc- 
tion can be identified. As an example, we illustrate in 
Fig. [2] the emerging asymptotic, averaged voltages across 
the first and the second junction for two different dimen- 
sionless temperatures D. We show how the parameter 
plane {ai, uo} is divided into regions of normal (i>i(ii) > 
and ^2(^1) > for /1 > ) and anomalous (^2(^1) < 
for Ii > 0) transport. In doing so here we do not dis- 
criminate between ANM and NNM. Both these transport 
behaviors are jointly presented. A more careful investi- 
gation of the parameter regimes, where negative mobil- 
ity appears, reveals a rich structure which is gradually 
smeared out by increasing temperature. 
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FIG. 2: The transport properties of the driven system of 
two Josephson junctions in the parameter space {ai,cj} of 
the ac driving acting on the first junction at a representa- 
tive coupling value of a — 0.77, the dc bias I\ — 0.05 and 
d2 = h = 0. In panels (a) and (b) the averaged voltage v\ 
across the first junction is shown at two dimensionless tem- 
peratures D = 5 • 10 -5 and 2.5 • 10 -4 , respectively. Panels 
(c) and (d): the averaged voltage V2 across the second junc- 
tion is depicted at the same temperatures as in panels (a) and 

(b) , respectively. Panel (e) presents an enlarged part of panel 

(c) : the plot reveals the island structure of regions of negative 
resistance of the second junction. 



We will now discuss the dependence of the voltage on 
the second junction of the dc bias I\ > acting on the 
first junction, in the regime of anomalous transport illus- 
trated in Fig. |2j We fix both the amplitude a\ and the 
frequency lj of the ac driving acting on the first junction. 
The result is presented in Fig. [3J Most prominent are 
the intervals of I\ where the voltage is negative. Two 
anomalous effects are detected: ANR for 1\ — » (the 
linear response regime) and NNR when 1\ has a value 
remote from zero. In the case presented in Fig.[3^a), one 
can observe three intervals of the bias I\ where NNR oc- 
curs. This is to be contrasted with the averaged voltage 
of the first junction which is never negative. However, 
one can note a rough synchronization in the dependency 
of v\ and V2 on the bias 1\\ simultaneous increases and 
decreases of both voltages in some intervals of the bias. 
Similar synchronization is also observed in the tempera- 
ture dependence, see Fig. |3j panels (b) and (c). A closer 
inspection of the panel (c) of Fig. [3] reveals another in- 
teresting result: there are two fundamentally different 
mechanisms generating anomalous transport. In the case 
I\ = 0.025, the negative resistance is induced by ther- 
mal fluctuations. There is a finite interval of tempera- 
ture where this effect occurs. For very low temperature 
(lower than 1.2 • 10 -5 ) and high temperature (greater 



0.20 
0.15 
0.10 
0.05 
0.00 
-0.05 
-0.10 

0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 



Vl,D: 

v 2 ,D- 



5 • 10" 5 
: 5 • 10" 5 



(a) 



0.015 



0.030 



0.045 



0.060 



._ I 1 = 0.025 
.... i x = 0.0456 




















s 





i " i 


" 6 10" 5 10" 4 10" 3 10 
D 




-0.06 
10 



FIG. 3: The long-time averaged voltages v\ and V2 of the 
first and second junctions, respectively. In panel (a) the de- 
pendence on the dc bias is depicted at a fixed temperature 
D — 5 • 10 -5 . Panels (b) and (c) depict the temperature de- 
pendence of V i and V2, respectively. In panel (c), please note 
two distinct mechanisms for the negative voltage of the sec- 
ond junction: noise-induced for I\ — 0.025 and deterministic, 
chaotic-assisted for ii = 0.456. The parameters are: the cou- 
pling strength a = 0.77, amplitude a\ — 1.7754 and frequency 
uj — 0.1876 of the ac driving, a2 = I2 = 0.0. 



than 2 • 10 -3 ) the voltage takes positive values. Be- 
tween these two temperatures, the voltage is negative. In 
contrast, in the case I\ = 0.0456, the negative resistance 
is generated by purely deterministic dynamics. Even for 
zero-temperature D = 0, the resistance is negative. For 
this chaos- assisted mechanism, temperature plays a de- 
structive role: increasing temperature monotonically di- 
minishes the negative voltage V2, and after crossing zero 
at some critical temperature, the voltage assumes posi- 
tive values. 

It should be emphasized that the anomalous transport 
effects discussed above are all caused by the coupling 
between two junctions - without coupling, the negative 
voltage vanishes. Fig. [4] shows how the average voltage 
of the second junction depends on the coupling constant 
a. We can see that there are always finite windows of a 
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FIG. 4: The influence of the coupling parameter a on trans- 
port properties of the second junction. The long-time aver- 
aged voltage V2 across the junction is depicted for two fixed 
values of the externally applied dc current I\ and D — 5-10 -5 . 
For I\ — 0.025 there are three intervals of a where the voltage 
is negative, while for ii = 0.0456 there are two intervals of 
such a. Other parameters are the same as in Fig. [3] 

for which this effect can be observed. The location and 
size of these windows depend on values of other system 
parameters. For the regimes depicted in Fig. [3J we illus- 
trate it in Fig. |4j In the case when the bias I\ = 0.025, 
up to coupling strength around a = 0.73 the average 
voltage is zero valued (numerically speaking the current 
is negligibly small). Starting from the mentioned value 
the transport settles in and later the voltage reaches the 
highest value of V2 ~ 0.18 at a « 0.741. Next, the voltage 
stays positive up to a ~ 0.753 and then changes its sign 
for the first time, resulting in an anomalous response - 
the voltage assumes negative values. This situation holds 
up to a « 0.758 and then the current turns into the posi- 
tive but very small values. The same scenario is repeated 
two more times - between the values of a ~ (0.76, 0.788) 
and (0.806,0.826). In the case h = 0.0456, there are two 
distinct intervals of the coupling constant for which the 
voltage V2 is negative. 

Adding a constant bias to the second particle has a 
destructive impact on the negative mobility - the higher 
the value of ^2, the smaller the parameter area where neg- 
ative mobility can be observed. For instance, we find 0.1 
to be the limiting value of I2 at which areas of negative 
mobility cease to exist everywhere in the analyzed area 
when 1\ = 0.05. Increasing I2 has a general smoothing 
effect on the transport properties of the second particle 
- at higher values of ^2, the fine details visible for I2 = 
gradually disappear. Increasing I 2 also causes all features 
to move slightly towards higher values of a\ and uo. It is 
therefore possible to find areas of the parameter space 
where I2 induces a negative mobility effect as well as ar- 
eas where it increases the strength of an existing negative 
mobility effect. 



It is of great importance to construct simple models 
which allow to explain and clarify the understanding of 
unusual and anomalous transport properties not only in 
physical but also in biological systems such as e.g. the 
bi-directionality of the net cargo transport inside living 
cells [7j. Such models can also serve as a basis for the 
construction of more realistic and quantitative models 
for transport of interacting carriers in collective systems. 
In this paper, we constructed a simple model of two cou- 
pled elements, described in terms of overdamped dynam- 
ics, which exhibits negative mobility both in the linear 
and non-linear regimes. There are regimes where both 
particles can be transported in the same direction as the 
external static force acting on the first particle. There are 
also regimes where only the second particle can be trans- 
ported in the opposite direction to bias. The interaction 
between two particles constitutes a crucial ingredient of 
the model: without coupling, the second particle is not 
transported at all, while the first particle can be trans- 
ported only in the direction of the dc bias. The system 
can be manipulated in other ways. The source of en- 
ergy, the driving I\(t) = l\ + cl\ cos(co>£), can be replaced 
by another source of energy such as e.g. an unbiased 
multi-harmonic force or, in biological systems, chemical 
reactions. An open question is what minimal conditions 
should be satisfied to generate anomalous transport, in 
particular bidirectional motion of two particles. Our con- 
jecture is that chaotic dynamics is necessary, but not suf- 
ficient. Finally, it is a promising topic which can stim- 
ulate experimentalists to perform measurements testing 
our findings in systems of two coupled Josephson junc- 
tions where coupling can be precisely controlled by the 
external resistance and, moreover, to plan experiments 
for other systems. 
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